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Abstract 


Multilayer designs of metasurfaces are required for practical 
applications because they have more capability to manipulate 
electromagnetic wave than single-layer designs. Here thus 
extend the Surface Susceptibility Model (SSM) integrated in 
FDTD from single layer to multilayers, consequently 
accelerating metasurfaces simulations with good accuracy. The 
SSM is integrated in FDTD by using Generalized Sheet 
Transition Conditions (GSTCs). To demonstrate the validity of 
the proposed algorithm, a numerical comparison is presented 
between simulation of two layer geometrical structures in 
conventional FDTD and two layer surface susceptibility model 
in FDTD-SSM. 
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Hundreds of thousands of 
subwavelength elements 


fabrication 


Extremely dense meshes and huge number of unknowns leads to 
enormous computational problems. 
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Modeling 


Inserting the model into 
numerical methods 
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Artificial electromagnetic surface 


Surface Susceptibility Model 


Material in nature 
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1. The unknowns are reduced from 9 to 3 for bi-axial case. 


2. Bi-axial case means that the structure do not generate cross polarization 
and transformation between electric to magnetic fields. 
14 


-1H$11- 


Surface susceptibility for bi-axial case 





The relationship between surface susceptibility and reflection coefficients: 
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Two sets of simulation is required to extract the surface susceptibilities for 
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Surface currents in FDTD 
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Bi-axial SSM in FDID 
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Flowchart of FDTD-SSM 
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update E in the whole computational space with conventional updating equations 
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Remarks: 


“Reduction of Meshing Lattice Size by using FDTD-SSM 


— FDTD 
4 FDTD-SSM:Al = 0.6 mm 
o FDTD-SSM:4A/ = 1.5 mm 
—-= FEDTD-SSM: A/ = 3 mm 





9.2 9.4 9.6 9.8 10 10.2 104 106 10.8 


freq (Ghz) 


1. The meshing lattice size of conventional FDTD is — —0.05 mm. 


2. As the meshing lattice size of FDID-SSM decreasing to 4 — 0.6 


50 


mm= 12 x 0.05 mm, its results overlap with that of conventional FDTD. 
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Remarks: 
1. 





A SSM simulation in FDTD-SSM is validated by a brute force simulation 
in HFSS. 


. If the distance between two layers is too small, the FDTD-SSM is not 


valid any more. 


. The FDTD-SSM is valid when the distance is bigger than one period. 
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A FDTD-SSM:Al = 0.6 mm 
© FDTD-SSM:A/ =1.5 mm 
--— FDTD-SSM:A/ = 3 mm 





9.2 9.4 9.6 9.8 10 10.2 104 10.6 10.6 


fre (GHz) 
TABLE I 
COMPUTATIONAL RESOURCES CONSUMPTION 
cell size time consum. 
Conv. FDTD | dx = dy = dz = 0.05 mm 2500 s 
FDTD-SSM dz = dy = dz = 0.6 mm IO s 
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TABLE I 
COMPUTATIONAL RESOURCES CONSUMPTION 


cell size time consum. 


Conv. FDTD | dz = dy = dz = 0.05 mm 38241 s 


FDTD-SSM | dz = dy = dz = 0.6 mm 12 s 
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Conclusion 





* A novel FDTD-SSM algorithm for accelerating metasurface simulation. 


e Surface susceptibility model, extraction, surface currents in FDTD 
* Advantages: simple & efficient 


* Applications of the FDTD-SSM algorithm 


* Reduction of meshing lattice size, effectiveness for different environment 


* Practical finite array, A six-layer structure 
* Consuming time 1s reduced from hours to seconds. 
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